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REVISION 3 CHANGES 


This document is a revised version of the previous one dated December 
1972 and supersedes that document. Major changes have been made in the guidance 
equations to increase the flexibility of the system and improve its performance 
characteristics . 

The guidance system now includes a transition mode following the normal 
entry mode. Significant changes have also been made in the constant heating rate 
mode and in the roll angle command logic. 


iii 



FOREWORD 


This document is one of a series of candidates for inclusion in 
a future revision of JSC -04217, "Space Shuttle Guidance, .Navigation 
and Control Design Equations". The enclosed has been prepared under 
NAS9-10268, Task No. 15-A, "GN & C Flight Equation Specification 
Support", and applies to functions 6 through 9 of the Entry Guidance 
Module (OG5) as defined in JSC- 03690, Rev. D, "Space Shuttle Orbiter 
Guidance, Navigation and Control Software Functional Requirements," 
dated January 1973. 
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INTRODUCTION 


The Entry-Guidance Routine presented here is designed to take the orbiter 
* vehicle from entry interface (h « 400,000 ft) through the critical heating phase of 
entry down to the start of the approach phase (h « 70,000 ft). The basic ideas are 
outlined in Refs. (1) - (5), (7). Simulation results demonstrating the feasibility of 
the concept are given in Refs. 2, 5, 6. 

There are four basic guidance modes: 

(1. ) An initial programmed- maneuver mode in which the 
vehicle is oriented with a constant roll angle and 
angle-of- attack. This mode is used until the specific 
force is greater than a pre-selected level. 

(2. ) A constant heating-rate mode during which the stag- 
nation-point heating rate is held constant at a pre- 
selected value, chosen essentially to minimize heat 
loads on the vehicle without violating maximum 
temperature constraints. 

(3. ) A reference trajectory mode during which the vehicle 
follows a pre- stored trajectory designed to get 
the vehicle to the terminal point with a minimum re- 
quired TPS weight, and without violating operational 
constraints on the vehicle. 

(4.) A transition mode during which the vehicle’s angle- 
of -attack is linearly ramped as a function of velocity 
from the back to the front side of the L/D curve. 

This mode also provides terminal conditions compatible 
with the start of an approach guidance system. 

Thermal control is provided by varying the magnitude of the roll angle so 
as to follow a density- vs. - speed profile. Density information is derived from IMU 
measurements of the aerodynamic specific force acting on the vehicle. A-priori 
knowledge of the vehicle’s mass, effective aerodynamic area, and drag coefficient 
<c D > are required in the process. 

During the constant heat rate and reference trajectory modes, range con- 
trol is provided by changing the angle -of -attack of the vehicle. Upper and lower 
limits on angle -of- attack are required in order not to violate operational constraints 
on the vehicle. Lateral trajectory control is obtained by reversing the direction 
of the roll angle. 
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In the transition-guidance mode, roll- angle magnitude is used to provide 
both altitude and range control. Limits on roll angle magnitude are used in order 
to alleviate potential problems due to sonic boom focal effects. 
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FUNCTIONAL FLOW DIAGRAM 


The basic information flow in the Entry Guidance Routine is shown in 
Figure 1. This is based on the guidance concept of Refs. 2, 5, 7. 

After the routine is entered, a series of targeting computations are made. 
This involves the computation of quantities such as the vehicle's current heading 
error (A lb ), and range to the target point (0 ) . 

The particular guidance mode to be entered is next determined. There 
are four possible guidance modes: 

(1. ) Initial programmed maneuver 

(2. ) Constant stagnation- point heating rate guidance 

(3. ) Stored reference-trajectory guidance 

(4. ) Transition guidance 


The constant heating-rate mode is entered when the measured specific force on the 
vehicle is greater than a preselected value. The reference-trajectory mode is 
entered when the magnitude of the vehicle’s relative velocity is less than a pre- 
selected value. The transition mode is also started on a relative -velocity criterion. 

In the programmed maneuver mode the vehicle is oriented with a constant 
roll angle and an angle-of-attack. This orientation is maintained until the heating- 
rate mode is entered. 

In the constant heating-rate mode, the density altitude (h D > required to 
attain the desired stagnation point heating rate (q^) is first computed. An angle- 
of-attack command (a ^ ) is next computed, based on the desired range-to-go for 
the heating- rate mode. Finally, roll- angle magnitude commands ( <j> ^ ) are com- 
puted to control the vertical- plane motion of the vehicle so as to follow a density- 
altitude vs. speed profile. No roll reversals take place in this mode. 

In the reference-trajectory mode, the required reference-trajectory 
quantities are first obtained from the stored table at the current speed. These in- 
clude angle-of-attack (« D ) , altitude (hp), range-to-go (r GD ), and the nominal 
heading error The angle-of-attack command (a^,) is then computed as a per- 

turbation from the reference value toj-j) based on the difference between the stored 
and measured values of range-to-go. Roll-angle magnitude commands are com- 
puted in the same manner as for the constant heating-rate mode, excpet that the 
desired density altitude (hj^) is from the stored table. Roll-angle direction is based 
on a comparison between the current heading error and the reference-trajectory 
value. 
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During the transition mode, desired values of angle- of- attack, density 
altitude, and allowable heading error are obtained from a stored trajectory as in 
the reference-trajectory mode. Range-to-go (r^^ ) is computed from a parabolic 
fit,relating nominal range to speed. The commanded angle -of- attack is pre- 
programmed as a function of speed, and during this mode a c is set equal to - 

Roll- angle magnitude commands are again computed in the same manner 

as for the constant heat rate mode except that the desired density altitude (h^) from 

the table lookup is modified by a term proportional to the vehicle T s range error 

(r*_ - r~ )• Roll angle magnitude and direction are also based on the vehicle 1 s 
GD Cj 

current heading error. If rj is not too large, roll reversals are performed as in 
the reference trajectory mode. If, however, t) becomes excessive, both the roll- 
angle magnitude and roll direction are commanded to reduce this error. 
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3 . 


INPUT AND OUTPUT VARIABLES 


m sm-e 


r 



Av 


At 

*p 


a c 


Input Variables 

Transformation matrix from stable-member to earth-fixed coordinates 
Vehicle position (stable-member coordinates) 

Vehicle velocity w. r.t. air mass (stable-member coordinates) 

IMU-measurement velocity change 

Time interval over which Av is taken 

Average value of angle-of-attack over previous At 

Average value of roll- angle over previous At 

Output Variables 
Angle-of-attack command 

Roll angle command 
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4. 

4.1 


DESCRIPTION OF EQUATIONS 


Targeting Calculations 
The range from the present location of the vehicle to the desired target 
point is computed from 


0 - cos U RE * J:rde ^ 


(1) 


where i„_ and i OTVir , are unit vectors at the present and target- point locations 

on the earth. 

The vehicle’s current heading error (A 0), defined as the difference between 
the current vehicle heading and the great- circle heading to the target point^is com- 
puted as: 

All) = cos 1 Qqq E • VRHE ^ Slgn ^—GCE X -VRHE* ' — RE^ ^ 

where v R jj E is the horizontal component of vehicle's relative velocity, and J.yRHE 
is a unit vector along ^rrr • The term ts a unit vector which is horizontal 

at the current vehicle location, and lies in the desired great-circle trajectory plane. 

4. 2 Initial Programmed- Maneuver Mode 

The roll-angle command (0 C > w.r.t. the relative velocity vector is for a 
constant roll attitude, i. e. 


0 C 0 CONST 


(3) 


The angle-of- attack command (« c ) is also maintained constant during this mode, 
i, e. 


a c "const 


(4) 


4, 3 Constant Heating-Rate Mode 

4.3.1 Thermal Control 

The required density <P D ) corresponding to a constant stagnation- point 
heating rate (q D ) is computed from the expression: 


P D ' C 1 q D / V RE 


6. 3 


(5) 
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where 

stant. 

tion: 


v RE is the velocity of the vehicle w. r.t. the air mass, and c^ is a con- 
This relation can be obtained from the stagnation- point heating- rate equa- 


0.5 

q = c 0 P 


3. 15 


V RE 


( 6 ) 


where q is the heating rate, p is the density of the atmosphere, and c Q is a 
constant. 

Thermal control is obtained by causing the vehicle to follow the desired 
density (p D > vs. speed (v RE ) profile. To accomplish this objective, an estimate 
of density at the current vehicle location is obtained from IMU-derived drag mea- 
surements. The basic relation for the density estimate (£) is: 


A 

P = 


2d 




(7) 


where m is the vehicle's mass, a is the effective aerodynamic area, and c D is 
the computed drag coefficient. The quantity v RE is the magnitude of the vehicle's 
velocity w. r. t. the air mass. The measured drag (d) is computed from IMU ve- 
locity-change date (Ay) using the relation: 


d = I Av • i VRE /At | 


( 8 ) 


where _iy RE is a unit vector along v RE ,and At is the time interval over which 
Av is measured. 


Rather than working with P D and p directly, it has been found to be 
more convenient to transform these quantities to equivalent density altitudes above 
the earth, i. e. h D and h ^ . In early simulations a simple exponential model was 
used to implement the transformation, based on the relation 


P = 



(9) 


where p„ is the sea-level density of the atmosphere, and h is its scale height. 
Currently a stored density- vs. altitude table for the 1962 standard atmosphere is 
used here. 


Using h„ and h as the primary guidance variables, a simple control law 

D P 

has been formulated which permits accurate control of density during the entire 
entry phase. The basic command is the desired change in vertical acceleration 

(Aa ) , which is computed from the relation: 

V D 


A a. 


v D " k 0l (h D _h p ) + k 02 (h D' h p )+ (h D" h l ) 


( 10 ) 
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* » « t 

where h^, , h D 


are the first and second derivatives of h^ and h^ . 


The quantities k^and are maintained as constant sensitivities 

during the constant-heat-rate and reference-trajectory modes. Prior to the start 
of the transition phase they are gradually increased to higher values which are 
then held constant during the transition mode. To obtain best performance in the 
presence of measurement noise, smoothing of the altitude-rate term (h^) and the 
nominal vertical acceleration (hp ) is required. 


The quantity hj is computed by back differencing the vehicle's altitude 
rate (hj). This altitude rate is computed using navigated information and a correc- 
tion term to account for the change in earth radius with changing geocentric latitude 
(Fischer -ellipsoid earth model). The quantity h^ is evaluated according to the 
following equations: 


*1 


= —RE * —RE ' r f 


(ID 


where r f is the time rate of change of earth radius due to the vehicle's latitude 

<P 1 

rate. r f is computed as: 


r pole ee cos( V sinU A ) *A 

[1-ee cos 2 ( l A ) ] l * & 


( 12 ) 


where r is the value of earth radius at the pole, ee is the square of the ec- 

c pole 

centricity of the Fischer-ellipsoid; and the terms & ^ are the vehicle s 

latitude and its time rate of change, respectively. 


The roll-angle command ( <?> c ) required to obtain Aa^ is obtained by 
commanding a change from the current roll angle ( dp), according to the relation 


(Aa v )/ | tan <p p \ 


0 c = *p- 


D 


(13) 


where the current vertical specific force (s^) obtained from the inertial mea 
surement unit. Special logic is provided to circumvent problems with Eq. (13) 
when 0 is very small. 
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4.3.2 Range Control 

Range control during the constant heating-rate mode is obtained by varying 
the angle- of- attack of the vehicle. The desired great- circle range to be covered in 
the constant heating-rate mode (r G. > is the difference between the present great- 
circle range to the terminal point ^ (r^ ) and the nominal range to be travelled 

in the final reference-trajectory mode (r r ), and in the approach- guidance rou- 

, . U REF 

tine (r^ ) , 1 . e. 

°APP 



- r. 


REF 


APP 


(14) 


It should be noted that the constant heating- rate mode is terminated at a selected 
final velocity (v p _ ) . 

q 

If various reasonable assumptions are made for the constant heating-rate 
mode, such as small flight-path angle, small changes in heading, small changes in 
altitude compared to the radius of the earth, and constant c D , then a simple closed- 
form solution is obtained for r Q . . The mathematical relation is: 



, 6. 3 rT 6. 3 , 

C 2 V RE " V F • ) 

q 


q c n (a/m) 


(15) 


where c 0 is a constant. By rearranging Eq. (15), the desired value of 

p » 

for the present q and required r ■ is obtained as: 

4 



y 6.3 6. 3 , 

c 2 V RE V F ■ -> 


q 2 r G , (a/m) 



) 


(16) 


The angle of attack command (a c ) required to achieve c D . is then com- 
puted as: ^ 


or. 


= Ot 


+ (c T 


- c„) k , 
D a' 


(17) 


where a p is the average value of angle-of-attack over the previous time step and 
k is a constant gain approximately equal to 3 c D / dor . 
a q 
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4. 3. 3 Lateral Control 


Roll reversals are not commanded during this guidance mode. The direc 
tion of the roll angle is chosen to reduce the vehicle’s initial heading error. 


4. 4 Reference Trajectory Mode 


4.4.1 Stored- Table Lookup 

The guidance concept makes use of a stored reference trajectory designed 
to minimize TPS weight for the maximum cross- range case of interest. The stored 
variables are angle -of- attack (<* N ), altitude (h N ), great-circle range-to-go to the 
target point (r GN ), and allowable heading error (tj n ). The independent or index 
variable is vehicle speed w. r.t. the air mass (y N ). In current simulation studies 
a 30-element stored array is used. 


A linear interpolation is used between the tabulated variables. For a 

typical case where the vehicle’s speed (v RE > lies between the stored values V N 

and v, T , a dummy variable £ is computed as: 

K+l 


K 


t ■ <’RE - V /( X + 1 ' V 


(18) 


The desired or reference- trajectory angle -of -attack (« D ) is then computed from: 


a D = (1 -»N + £ a N 
u 1 k+1 


( 10 ) 


Similar relations are used to find h Q , r^^ , and r) D . 

4. 4. 2 Thermal Control 

The basic idea is to roll the vehicle about its velocity vector so as to 
follow a stored density vs. speed profile. The roll- angle magnitude computations 
are essentially the same as for the constant heating -rate mode, except that the 
desired density altitude (h D ) is obtained from the stored reference table (rather 
than from Eq. 9). 

4. 4. 3 Range Control 

Angle -of -attack commands ( a c > are computed as a combination of the 
desired and reference-trajectory range-to-go (r G and r GD ). The basic relation is 

“C = a D + k ff (r G ' r GD^ v RE (20) 
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The multiplying factor associated with the range -difference term is varied inversely 
with vehicle velocity (k^/v^p ). The maximum value of angle -of- attack is limited 
to a preselected level . 

4.4.4 Lateral Control 

The lateral control of the vehicle f s trajectory is accomplished by roll 
modulation, i. e. reversing the direction of roll. For this purpose it is necessary 
first to see if the present roll angle ( 0 ) is directed so as to increase the 

heading error (A l/> ). This is accomplished by looking at the test quantity de- 

fined by: 

% = A ^OLD (21) 


where 0 q LD is the preceding roll-angle command. If it turns out that is 

tending to increase A 0 ,then a check is made to see if the magnitude of the current 
heading error ( rj ) is greater than a specified fraction of the reference-tra- 




jectory value (rj^) . If it is then seen that TJ is greater than k 
the roll- command direction is reversed from the previous command 

4. 5 Transition Mode 


< 


OLD 


then 

). 


4.5.1 Nominal Trajectory Description 

The transition mode uses a stored trajectory as does the reference-trajec- 
tory mode. The nominal trajectory for this phase is not based on temperature or 
TPS considerations, but rather on providing maximum footprint capability (to null 
post-blackout navigation errors) while avoiding violation of vehicle constraints 
(eleven deflections, hinge moments, g loads, etc.) 

Angle -of- attack, altitude, and allowable heading error are stored as a 
function of speed as in the reference-trajectory mode. Nominal range to the target 
^ r GN^ h° wever > computed from a parabolic function of range vs. speed, i. e. 

r GN = ^ (v RE “ V F* + c R 2 <v RE “ V F* (22) 


where v-p is the value of velocity at which the guidance mode is terminated. This 
form is used for the nominal range computation rather than a table lookup because 
it provides smoother control histories. 
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4.5.2 Angle-of- Attack 

Angle -of -attack is pre-programmed as a function of velocity during transi- 
tion. This technique accomplishes the three following objectives; (1) keeps the 
vehicle in a controllable orientation ( a vs. Mach no. ), (2) helps avoid violation of 

elevon deflection and hinge-moment constraints, and (3) keeps the vehicle flying at 
near maximum L/D (thus providing good maneuver capability). 

4.5.3 Range Control 

Range is controlled by altering the altitude -vs. -velocity profile flown by 
the vehicle. An increase in altitude provides a lower density (or drag) thereby 
increasing range. The opposite is true for a decrease in altitude. The altitude-vs.- 
velocity level at which the vehicle is commanded to fly is computed from the 
nominal reference altitude and the vehicle's current range error as follows: 


h D = h N +k(v) (r G ‘ r GN } 


(23) 


where k(v) is a gain that is linearly increased with decreasing velocity. This 
value of desired altitude is then used to compute roll- angle magnitude commands 
as in the constant heat rate and reference-trajectory modes. During this mode 
the limits on maximum roll-angle are modified to keep the vehicle from violating 
sonic boom over-pressure constraints. 


4. 5. 4 Lateral Control 


Lateral control during transition is accomplished by use of a double- dead- 
band logic. These deadbands are defined by the quantities k and k (Note: 

V l V 2 

k > k ). If the vehicle's current heading error rj is less than k r? 

W 2 1 or TJ ' 

the roll reversal logic of the reference- trajectory mode is used for lateral 
control. 


D ' 


If the vehicle's heading error is greater than k^ rf^ } and if the vehicle 
is currently too far from the target point (i. e. , ) * then both roll- angle 

magnitude and roll- angle direction are commanded to decrease T) . This is accomp- 
lished by commanding a roll- angle magnitude of 45° (if the normally commanded 
roll-angle is less than 45°) with the proper sign to reduce the heading error. The 
final output roll command is first limited by the maximum roll-angle constraints. 

A range-error criterion is used in this lateral logic because in those 
cases where the vehicle is too far from the target, the normal roll angle commands 
are small (in order to increase in-plane L/D to increase range). These small 
roll angles tend to decrease the vehicle’s lateral capability, and the vehicle would fly 
past the target if some explicit lateral control were not employed. 
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5* DETAILED FLOW DIAGRAMS 

This section contains detailed flow diagrams of the Entry Guidance Routine. 

program constants 



Figure 2a. Entry Guidance Routine, Detailed Flow Diagram 


5-1 







Figure 2b. Entry Guidance Routine, Detailed Flow Diagram 
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Figure 2c. Entry Guidance Routine, Detailed Flow Diagram 
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OUTPUT VARIABLES, 



Figure 2d. Entry Guidance Routine, Detailed Flow Diagram 


5-4 






Figure 2e. Entry Guidance Routine, Detailed Flow Diagram 
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Figure 2f . Entry Guidance Routine, Detailed Flow Diagram 


5-0 











Figure 2g. Entry Guidance Routine, Detailed Flow Diagram 
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Figure 2h. Entry Guidance Routine, Detailed Flow Diagram 


5-8 









Figure 2i. Entry Guidance Routine, Detailed Flow Diagram 
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(Fig. 21) 


Figure 2j. Entry Guidance Routine, Detailed Flow Diagram 
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Figure 21. Entry Guidance Routine, Detailed Flow Diagram 
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